Magnetic field of Josephson vortices outside superconductors 



00 
O 

o 

(N 



V. G. Kogan 1 and V. V. Dobrovitski 1 

1 Ames Laboratory and Department of Physics and Astronomy, ISU, Ames, I A 50011 

(Dated: March 20, 2008) 

We consider the structure of Josephson vortices approaching the junction boundary with vac- 
uum in large area Josephson junctions with the Josephson length Xj large relative to the London 
penetration depth Xl- Using the stability argument for one-dimentional solitons with respect to 
2D perturbations, it is shown that on the scale Aj the Josephson vortices do not spread near the 
boundary in the direction of the junction. The field distribution in vacuum due to the Josephson 
vortex is evaluated, the information needed for the Scanning SQUID Microscopy. 

PACS numbers: 74.50.+r, 74.78.Na 
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A new experimental technique, the scanning SQUID 
microscopy (SSM), has recently been developed for mea- 
suring magnetic field distributions due to vortices exit- 
ing superconducting samples. 1 Knowing the distributions 
one can, in principle, extract the London penetration 
depth Xl (either isotropic or anisotropic) and its tem- 
perature dependence. 2 

To extract the superconducting parameters from the 
measured field above the sample one has to solve the 
problem of the field created by a vortex terminating at 
the sample surface inside and outside the sample. For 
isotropic bulk materials the solution has been given by 
J. Pearl who utilized the cylindrical symmetry of a vortex 
normal to the sample surfaced For the anisotropic case, a 
more general approach should be employed to match so- 
lutions of London equations inside and of Maxwell equa- 
tions outside^ The most likely situation of a vortex per- 
pendicular to the sample surface was described and im- 
plemented in analysis of SSM data in Ref.[f| These cal- 
culations showed that vortices spread out in the super- 
conductor while approaching the surface. For the vortex 
axis along z, the transverse components h x ,y appear de- 
scribing field lines bending out of the vortex axis in the 
subsurface layer of a depth on the order of Xl ■ 

Originally, the SSM technique has been designed for 
studies of Josephson boundaries between misoriented 
high-T c superconducting crystallites. Some of the most 
convincing demonstrations of the d-wave symmetry of 
these materials were obtained by using SSM to detect 
spontaneous half-flux-quantum vortices at the intersec- 
tion of three Josephson boundaries between crystallines 
of a proper misorientationi^ Other configurations in- 
volving Josephson boundaries have been shown to con- 
tain spontaneous magnetic flux (e.g., closed triangular 
and hexagonal boundaries^ faceted boundaries with al- 
ternating critical current^). However, absence of theo- 
retical description of Josephson vortices terminating on 
the sample surface prevented a reliable interpretation of 
a large body of SSM data. 

The Josephson length Xj for the boundaries in ques- 
tion (the scale of the field distribution in the Josephson 
vortex that plays a role of Xl for the bulk vortices) might 
be on the order of microns. This suggests that a consid- 



erable spreading effect may take place at the distance of 
order of Aj near the surface. Such an effect, if exists, 
should be taken into account when relating the outside 
field distributions to the internal structure of Josephson 
vortices (J- vortices). 

We argue in this communication that Josephson vor- 
tices do not spread in the direction of the junction in the 
limit 



X L <C Xj 
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This conclusion makes the SSM data analysis not only 
possible but quite straightforward. 

Let us start with an infinite two-dimensional (2D) 
Josephson junction in the plane (a;, z) between two su- 
perconducting banks occupying the half-spaces y > 
and y < 0. The stationary J-vortex directed along z in 
such a junction is described by the sine-Gordon equation 
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where f{x) is the gauge-invariant phase difference ] 10 ' 11 : 12 

— 1/2 

The Josephson length Xj oc j c may vary due to dif- 
ferences in crystals misalignment and in quality of the 
boundary (J c is the critical Josephson current density). 
The primes denote differentiation with respect to x. One 
follows the standard procedure: multiply (J2J) by <p' to get 
the first integral X 2 j<p n /2 = C — cos if. The Josephson 
current and the field (oc if') at the junction should van- 
ish as \x\ — » oo; this yields C = 1 for 93(00) = 2tt. We 
then obtain: 



4tan _1 (e a/Aj )- 



The field at the junction plane: 



h z (x) 
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0o = irhc/\e\ is the flux quantum (the thickness of the 
insulating layer is assumed small relative to Xl)- The 
vortex described by Eqs. (|3|) and ((4|) is infinite in the 
z direction; the problem and the solution are, in fact, 
one-dimensional. 

Let us now remove the superconductors from the half- 
space z > 0; thus the plane z = is a free boundary 
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with vacuum. The question arises what changes the J- 
vortex (occupying now only z < 0) should undergo? The 
problem is no longer uniform in the z direction, and the 
z independence of the solution of Eq. cannot be as- 
sumed. In other words, we have now ip = p{x 1 z) and the 
sine-Gordon equation in two spacial dimensions: 

Aj VV = sin<^ , (5) 

V 2 is the 2D Laplacian. We look for solutions at the 
half- plane — oo < z < describing a J- vortex somehow 
spreading as z — > —0. The solution should satisfy certain 
boundary conditions: at z — > — oo, tp must be of the form 
([3]). Besides, at any z < 0, the vortex field should vanish 
far from the vortex "core" : 



dip 
dx 



= 0. 
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This formidable, at first sight, problem does have an 
exact and unique solution due to the "nonexistence" the- 
orem by Derric k 13 i 14 that states that there are no stable 
static solutions of the sine-Gordon equation in more than 
one spatial dimension. 

We reproduce this argument for the 2D case of interest 
here. The energy functional E{p(x, z)}, which generates 
Eq. (O, can be written in the form: 



E 



dx dz 
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K + U. 



(7) 



where we have dropped a constant prefactor irrelevant 
for the following. Let us assume that a function (p(x, z) 
exists such that 6E/5<p = 0; in other words, that the 
sine-Gordon equation has a 2D bounded solution. We 
now evaluate the energy for the phase ip(rjx, rjz) with an 
arbitrary scaling factor r\. Substitute this in the integral 
([7]) and change integration variables to r]x, r\z to obtain: 



E(r)) = K + U/rf 



(8) 



where the kinetic and potential parts, K and U, are 
defined for the original phase p(x,z). Since p(rjx,r]z) 
is a solution of 5E = for r\ = 1, we must 
have dE(ri)/dr]\ n=1 = 0. Equation ([8]) however gives 
dE{rj) / dr)\ v =i = —2U < 0. Thus, there are no 2D static 
solutions-cither stable or unstable-of Eq. ([5]) . 

Note that for ID case E(ri) = i]K + U/rj, 
dE(r))/dr]\ v= i = K — U and ID solutions with K = U do 
exist and are stable. 

It should be noted that the above argument implies 
an infinite plane x,z; both integrals ([7]) are taken on 
the interval (—00,00). For the junction of our interest 
— 00 < z < 0. Clearly, the scaling employed above can 
still be used for the half-plane provided the coordinate z 
is counted from the junction edge so that the edge is left 
in place under the scaling transform. 

The variational minimization of the energy ([JJ defined 
on the half-plane involves the integration by parts which 
yields the boundary condition 
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~d~z 



0, 



(9) 



which is invariant with repect to the scaling employed 
above. This condition translates to h x (x,0,0) — along 
with g z {x, +0, -0) = g y {x, +0, -0) = 0. 

We thus conclude that Eq. §5§ with z independent ip(x) 
is the only solution possible. In other words, the static 
Josephson vortex does not spread approaching the sam- 
ple surface. 

It is worth mentioning that dynamic 2D solutions of 
the time dependent sine-Gordon equation do exist. A 
few examples are given in Ref.[l5| where the initial es- 
sentially two-dimensional soliton-like structure has been 
shown (numerically) to evolve in time toward a straight 
ID soliton. 

We stress that the above conclusion holds when effects 
on the scale of London penetration depth are neglected. 
Of course, at distances of the order Xl from the surface, 
the field deviates from the z direction and some spreading 
takes place. This effect, however, is hardly relevant for 
the SSM method because the low bound of the spatial 
resolution of the SQUID loops is set by the loop size, 
which exceeds by much the London depthi 

Also, one has to keep in mind that near the free surface, 
the stray fields outside the superconductor may affect the 
very equation governing the behavior of the phase differ- 
ence p. This question has been discussed for thin-fihn 
junctionS ) 16 i 17 ' 18 where it was shown that the governing 
equation replacing the "bulk" Eq. ([5]) becomes integral, 
in other words, the stray fields make the problem nonlo- 
cal. This perturbation, however, penetrates the system 
only down to depths on the order of Xl and can be dis- 
regarded in our approximation. 
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FIG. 1: Contours of h z (x,y) = const for Xj = 1 and 20 = 0.5. 
The horizontal axis is x along the junction. The field in units 
of 4>q/A-k 2 X 2 j reaches maximum at the origin h z (0, 0) = 3.507. 
Values of h z are indicated on the contours. 



Since the field h in vacuum satisfies curl h = div h = 0, 
one can introduce in the upper half-space z > a scalar 
potential: 



vV = 0. 



(10) 
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In magnetostatic terms, this corresponds to a linear 
"charge" with the density 
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FIG. 2: Same as Fig. [T] for Zo = 1. The maximum field is 
h z (0,0) = 1.474. 
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Therefore, the potential ip reads: 



where 



dx 



J J-oo R cosh(x/Aj) ' 



R 2 = (x - x Q ) 2 + y 2 + z, 
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The field component h z , which is actually measured by 
SSM, is readily obtained: 



h z (x Q ,y Q ,z ) 



<t>0 zo 
47r 2 A/ 



dx 



R 3 cosh(x/Aj) 



(15) 
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FIG. 3: Same as Fig.[T] for zq = 2. The maximum field is 
M0,0) = 0.5388. 



This integral is well convergent at any finite z and is easy 
for numerical evaluation. It diverges at the junction line 
at z — 0; the divergence, however, should be truncated 
at distances of the order A^ which are disregarded in our 
model. 

Contours of h z (x,y) = const are shown in Figs. [HO 
For this calculation we set Aj = 1 and z$ = 0.5, 1 and 2. 

In actual SSM, the field h z (x,y;z) at a fixed height 
z should be integrated over the SQUID area. In some 
implementations it is more convenient to have the 2D 
Fourier transform h z (k x , k y ] z) instead of the real space 
distribution. After straightforward algebra we obtain: 



h z {k x , ky] z) 



-kz 



cosh(7rfc x Aj/2) 



. k=Jkl + kl. (16) 



The boundary conditions for the Laplace equation are 
provided by the Josephson field "sticking out" of the 
plane z — from a belt of the width 2A^ along the junc- 
tion. With the help of Eq. (U|) we obtain: 



dz 



= 2A L %)— p-sech-. (11) 



In particular this shows that the total flux through any 
plane z = zq is h z (\t = 0; zq) — 4>q. 
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